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ANALYTICAL SOLUTION OF SPACE-TIME FRACTIONAL 
TELEGRAPH-TYPE EQUATIONS INVOLVING HILFER AND 
HADAMARD DERIVATIVES 

RAM K. SAXENA ROBERTO GARRA^, AND ENZO ORSINGHER^ 


Abstract. In this paper we consider space-time fractional telegraph equations, 
where the time derivatives are intended in the sense of Hilfer and Hadamard 
while the space fractional derivatives are meant in the sense of Riesz-Feller. We 
provide the Fourier transforms of the solutions of some Cauchy problems for these 
fractional equations. Probabilistic interpretations of some specific cases are also 
provided. 


1. Introduction 


In recent years, an increasing interest for the analysis and applications of space and 
time-fractional generalizations of the telegraph-type equations has been developed 
in the literature. One of the first works in this direction was the paper by Orsingher 
and Zhao na about the space-fractional telegraph equation and then a number of 
papers regarding applications in probability mm, in physics ID El or merely of 
mathematical interest mu appeared in the literature. 

In the first part of this paper we study the fractional telegraph equation involving 
Hadamard-type time-fractional derivatives, that is 


( 1 . 1 ) 




u + 2X 



2d‘^u 


with x G M and t > to > 0. We denote with the symbol u the Caputo-type 

modihcation of the Hadamard derivative of order u recently introduced in [TU] as 
follows 


( 1 . 2 ) 






for n — 1 < 1 / < n and n G N. We will explain the reason why we choose to use this 
symbol for the Caputo-like Hadamard fractional operator and recall some of its main 
properties in the next section devoted to mathematical preliminaries. Within the 
framework of the analysis of the fractional telegraph equation involving Hadamard 
derivatives we will present analytical results concerning the Fourier transform of the 
fundamental solution and give some insights regarding the probabilistic meaning of 
the obtained results. 
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In the second part of the paper, we consider the generalized telegraph equation 
involving Hilfer fractional derivatives in time and Riesz-Feller fractional deriva¬ 
tives xDq in space. The generalized telegraph equation is the following partial 
differential equation 


(1.3) 


dt"^ dt 


— uu, 


a; > 0 


and some fractional generalizations of this equation have been recently studied in 
a number of papers (see for example EE] and the references therein). It is well- 
known that equation 01.31) arises, for example, in the analysis of transmission lines 
with losses. 

Here we give a new useful representation formula for the solution of the Cauchy 
problem for the equation 


(1.4) D‘P’^u{x, t) + 2\D'l'^u{x, t) = (? xDgu{x, t)—uu{x, t) + F(x, t), x G M, t > 0 


We observe that general results about non-homogeneous generalized fractional tele¬ 
graph equations involving Hilfer and Riesz-Feller operators have been discussed in 
im. Here we focus our attention on a specihc case, that of equation fll.41) . where the 
solution can be expressed in terms of linear combination of two-parameter Mittag- 
Leffler functions. This permits us to give an interesting probabilistic interpretation 
of the solutions when the parameters 7 and 6 are suitably chosen. We also show that 
the obtained results extend those appeared so far for the classical and the fractional 
telegraph equation. 


2. Preliminaries 

For the reader’s convenience we hrst present the necessary mathematical prelimi¬ 
naries about fractional operators and their main properties used in the next sections. 
We moreover give a brief recall about the classical telegraph process and its 

relationship with the telegraph equation. 


2.1. Fractional differential operators. The classical Hadamard fractional deriv¬ 
ative reads (see for example [16]) 


( 2 . 1 ) 






n—u—l 


dr 

u{t)—, 


where f > to > 0 and n — 1 < 1 / < n, with n G N. 

In the recent paper the Caputo-type regularized counterpart of the Hadamard 
fractional derivative has been introduced as 


( 2 . 2 ) 






where t > to > 0 and n — 1 < < n, with n G N. In this paper we adopt the symbol 

for the Caputo-type Hadamard fractional derivative in order to highlight the 
relationship with the differential operator 5 = t4. Indeed, we observe that, for z/ = 1 
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this operator coincides by definition with 6. One of the useful properties that we 
will be employed in the next section is the following one 

( 4 ) 4 ^) = r 4 +t 4 ) (' 4 ) ’ >-1/m 

which can be proved by simple calculations. We remark that, in the case /3 = 0, as 
can be simply understood in view of the definition (2.2), we have that 


(2.4) 


/ iiy 

t— const. = 0. 

V */ 


We refer to [IH] for a complete discussion about the functional setting and main 
properties of this integro-differential operator. 

In a series of works (see [9] and the references therein), Hilfer studied applications 
of a generalized fractional operator having the Riemann-Liouville and the Caputo 
derivatives as specific cases. The Hilfer’s derivative is defined as 




r(z/(n-/i)) 


dy (t - y) 






dy- 


r((l-i^)(n-/i) 


u{x, z)dz 


(2.5) 




U e [0,1), P e (0,1). 


Clearly n — 1 < (1 — z/)(?7, — /i) < n and is the Riemann-Liouville fractional 

derivative and is the Riemann-Liouville integral (see e.g. [HE]). Hereafter 

we will take for simplicity a = 0. 

The Laplace transform of the Hilfer fractional derivative w.r. to the time reads 
(see [T^ l 


(2.6) C(yD^+u) {x, s) = / e ‘^^D^+u{x,t)dt 


72—1 


s>^u{x, s)-J2 t) 


k=0 


dt^ 


t=o 


where (£u) (s) = u{x, s). 

We now recall the definition of the Riesz-Feller fractional derivative and for more 
details we refer, for example, to the encyclopedical book by Samko et ah [16]. For 
0 < a < 2 and \6\ < min(Q!, 2 — 0 ;), the Riesz-Feller derivative is defined as 

(2.7) 


UDSf) (x) = 

r(l-a)f . 7r(a + «) f" f(x + 0 - f(x) , . x(a - 0) r f(x - 0 - f(x) 

For 6 = 0, the Riesz-Feller derivative becomes the Riesz fractional derivative, that 

is 


d^ 



a/2 


( 2 . 8 ) 


(m;/) (a:) 
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The Fourier transform of fl2.7p reads 

(2.9) (/3) = -\l3re‘^’“^0rW) = f (fi), 

where = (Tf) (j3). 

We note that the fundamental solution to the space-fractional differential equation 
involving the Riesz-Feller derivative 

has Fourier transform 


( 2 . 10 ) 


p*{(3,t) = exp{ — \(3\°‘te'^'’^ sign(3}, 


which is the characteristic function of a stable process SaicTyj, p;t), with p = 0 , 
a = cos ^ and 7 = (see e.g. ISO])- 


2 . 2 . The telegraph process. The so-called telegraph process, T(f), f > 0, de¬ 
scribes the random motion of a particle moving on the real line with hnite velocity c 
and alternating two possible directions of motions (forward or backward) at Poisson 
paced times with constant rate A > 0. This simple finite-velocity random motion 
was firstly suggested by the description of particles transport. In this framework, 
the first derivation given in literature, as far as we kwow, was given by Fiirth in 
a discussion of several models of fluctuation phenomena in physics and biology. A 
similar model was also considered by G.I.Taylor [T2] to treat turbulent diffusion 
and then studied in detail by Sidney Goldstein in |H] and Mark Kac in m- This 
process is called in the literature telegraph process because the probability law of 
T{t) coincides with the fundamental solution of the hyperbolic telegraph equation. 
Moreover, as a limiting case it coincides with the transition density of the classical 
Brownian motion. 

In more detail, the classical symmetric telegraph process is dehned as 

( 2 . 11 ) T{t) = 1/(0) [ (-l)^(")ds, t > 0, 

Jo 

where V (0) is a two-valued random variable independent of the Poisson process 
Af{t), t > 0 . 

The component of the unconditional distribution of the telegraph process concen¬ 
trated inside the interval (—cf, +ct), is given by 

( 2 . 12 ) 

e-^ 

P{T{t) G dx} = dx—— 

The component of the unconditional distribution that pertains to the Poisson prob¬ 
ability of no changes of directions is concentrated on the boundary i.e. x = ±cf. 



A/n f 


X" 


d 


A 


+ T^Io - X‘ 

ot \ c 


|x| < ct. 


(2.13) 


P{r{t) = ±ct} 


2 
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Hence we are able to give in explicit form the density f{x,t) of the distribution of 
T(t), that is 


.-At r 




+ 


2c 

At 


A 


A/n -Vc2^ 




9 r 


-[5{ct — x) + 5{ct + x)], 



l(-ct,+ct) {x) 


where ©(■) is the Heaviside function and (5(-) is the Dirac delta function. We remark 
that the component of the distribution of the telegraph process concentrated in 
(—ct, +cf), given by formula fl2.12p . is the solution to the Cauchy problem 


(2.14) 


< 


|f + 2Af = c 

p(x,0) = (5(x), 


2 d^p 
dx^ ’ 


dp 

dt 


(x,t) 


= 0 . 


t=o 


3. The fractional telegraph equation involving Hadamard 

FRACTIONAL DERIVATIVES 

In this section we study the following time-fractional telegraph equation involving 
Caputo-like regularized Hadamard derivatives introduced in the previous section 

(4) (4) “+^44) 

It corresponds to a generalization of the telegraph equation with variable coefficients. 
The main results presented in this section are strongly inspired by the previous one 
obtained by Beghin and Orsingher in [T3]. Indeed we will do an Ansatz about the 
form of the Fourier transform of the fundamental solution that is based on the results 
proved in [13] and then we verify our assumption. We will also explain the heuristic 
reasonableness of our statement. 


Theorem 3.1. The Fourier transform of the solution to the Cauchy problem 

+ = z/e (0,1), to>0 

u(x,to) ^ (5(x), 


(3.2) 


is given by 


=0 

dt li=io ’ 


(3.3) 


u*{P,t) = 


2 

+ I 1- 


1 + 


A 


A 


Eu,i [piln 
Eu,i ( 


where 

(3.4) 


r]i = -\ + T]2 = -X- -x/A^ - c^l3 
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Proof. Let us introduce the operator 


n,= 




+ 2A 



By using the property fl2.3p of the Caputo-like Hadamard derivative, we have that 



E, 


i/,i 




T{uk+1) {ln{t/to)Y^ '"Vi 
V{i>k + l — p) V{i/k + l) 


n^k' + i) 


ViEu,i 



1 


where we used the fact that the Caputo-type Hadamard derivative of a constant is 
zero. On the basis of this result we have by straight calculation that 


(3.5) 




as claimed. 


□ 


Result (3.3) is heuristically motivated by the fact that the Fourier transform of 
the solution of the fractional telegraph equation involving Hadamard derivatives can 
be obtained from the solution of the fractional problem involving Caputo derivatives 
by means of the deterministic time-change f —)■ In f. 

This observation implies that, for u = 1/2, formula 03.31) becomes the charachteris- 
tic function of a telegraph process at a time-changed Brownian motion, as suggested 
by the results of section 4 of [13]. In fact, for u = 1/2, formula 03.31) can be written 
as 


(3.6) 


A 


poo _^ 




e ^ 


poo _2“ 


-A^ r 




dz 


—Xz 


e *0 


that coincides with the characteristic function of the process 


2-y/A2-c2/32 g-zy'A2-c2/32 


dz 


(3.7) 


w(t) = r 



1 


that is a telegraph proces with random Brownian time with a logarithmic determin¬ 
istic time-change (see [T3]L 

We moreover observe that for —)■ 1, we have studied a fractional generalization 

of the telegraph equation with variable coefficients 


(3.8) 


A 


. , d\ 2 d'^u 

« + 2A|(^ u = 


Then, it is simple to prove that the fundamental solution of this telegraph-type equa¬ 
tion coincides with the law of the classical telegraph process with the aforementioned 
deterministic time change. 
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On the basis of the previous considerations, we now consider the space-time frac¬ 
tional telegraph equation 


(3.9) 




involving time-fractional Hadamard derivatives and space-fractional Riesz derivative 
(see the previous section). Equation fId.lOp has been previously studied in |1] in the 
case where the Hadamard derivatives are replaced with Caputo fractional derivatives. 
In this paper the authors have shown that the fundamental solution of the space-time 
fractional telegraph-type equation 


(3.10) 


, d^u 

-h 2A- 

dt^'^ dt^ 


—c 


V dx^) 


e (0,^], a e (0,2] 


coincides with the distribution of the composition of a a stable process S^(t), t > 0 
with the positively-valued process 

(3.11) C'it) = inf{s > 0 ; ?{'^(s) = + {2XY''H^{s) >t}, t> 0, 

where and independent positively skewed stable processes of order 2u 

and u respectively. With the next proposition we generalize Theorem 4.1 of [1] to 
the case of equation fl3.10l) 

Proposition 3.2. For v G (0,1/2] and a G (0,2], the fundamental solution of 
equation 03.1011 coincides with the probability law of the process 


(3.12) 


W{t) = R" 1 ^ > 0 


io 


and has Fourier transform 


(3.13) 

where 

(3.14) 


u*W,t) = - 


+ 1 


1 + 






E,i ( ri/n" - 
. 2 \R\a I to 


Eu,i ( rs/n" ) , 


ri = -\ + 


= -A- 


2 -C2|d|« 


We neglect the complete proof of this proposition because it consists of a simple 
combination of the arguments of the proof of Theorem 4.1 in [1] and of the previous 
considerations about the role of Hadamard time-fractional derivatives that induces 
a deterministic logarithmic time-change. However we remark that this result can 
be easilty generalized to the multidimensional case (where the fractional Laplacian 
appears) and gives a more general probabilistic interpretation of the fundamental 
solution of the space-time fractional telegraph equation. 
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3.1. Decomposing the fractional telegraph equation involving Hadamard 
derivatives. Here we show that the fractional equation fl3.ll) can be easily found 
by decoupling it into a system of two fractional equations and we provide a brief 
discussion about the utility of this equivalent formulation. By simple calculation it 
can be proved that fl3.ip emerges from the system 


(3.15) 


{t§^y w + 2Xw =-c^ 


The utility of the decomposition fl3.15p of equation fl3.8p is twofold. It gives useful 
insights on the meaning of the equation fl3.ip both from the probabilistic and physical 
points of view. From the physical points of view this is a generalization of the 
fractional Cattaneo heat model widely studied in the literature (see for example [15] 
and the references therein). The generalization stems from in the fact that, replacing 
the Caputo derivative with the Hadamard derivative, we obtain a deterministic time- 
change leading to a slower, logarithmic in time, heat wave propagation. On the other 
hand, we recall that the probability law of the telegraph process can be obtained by 
solving the coupled system 


(3.16) 

(3.17) 

where 


d_ 

dt 

d_ 

dt 


dw 


u = 


dx 

w + 2Xw = 



(3.18) u{x,t) = + b{x,t), w{x,t) = — b{x,t), 

and f{x,t) represents the probability that the particle is near x at time t with 
forward velocity and b{x, t) the probability that the particle is near x at time t with 
backward velocity. 


4. The generalized fragtional telegraph equation involving Hilfer 

DERIVATIVES 

In this section we study the generalized non-homogeneous fractional telegraph 
equation involving Hilfer derivatives in time and Riesz-Feller derivatives in space 

(4.1) Dy’^u{x, t)+2XD^’^u{x, t) = (? xDqu{x, t)—uju{x, t)+F{x, t), a; G M, t > 0. 

In this framework our hrst general result is given by the following 

Theorem 4.1. Consider the generalized non-homogeneous space-time fractional tele¬ 
graph eguation fl4.ll) subject to the following constraints 

(4.2) 0<7<1, 0<(5<1, 0<q;<2 |0|< min(cK, 2 — 0;). 

If the eguation fl4.ll) is subject to the following initial and boundary coinditions 

(4.3) 0) = fi{x), 


(4.4) 


( 1 _ 5 )( 2 - 27 ) 

dt 


u{x, t) 


= 0 , 


t=o 
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(4.5) 

(4.6) 

then its solution reads 

(4.7) 

, , 1 
u[x, t) = 


'o+ 


m(x,0) = f2{x), 
lim u(x, t) = 0, 

|3:|—>-oo 


H-oo 


27r(^ - T]) 


- ^7,7+25-27^-1 (^^^)}] Fl d(3 


"+00 


27r(^ - r]) y_, 

^ 27r(^ -V) lo ' 

where 

(4.8) 

and b = u + c^il^a- 


^7-1 


[f< M^7,7+5-75- ^7,7+5-75/a (/3)e(rp(-i/3a:)d/3 

/ + 00 

F*{I3, t — t) exp{—ij3x)dj3dT^ 

-oo 

^ = -A + r] = -x- \/P^, 


Proof. If we apply the Laplace transform with respect to the time variable t and 
Fourier transform with respect to the space variable x, by using the initial conditions 
and formulas fl2.6p and fl2.9p . we obtain 

(4.9) s^'^u*{F s) - /*(/?) + 2A [s'^u*{F s) - 

= s) - uju*{/3, s) + F*{l3, s), 

where the symbol * represents the Fourier transform with respect to the space vari¬ 
able X and u(/3) stands for the Laplace transform with respect to the time variable. 
From fl4.9p . we have that 


(4.10) 


u*(/3,s) = 


s 1 -^( 2 - 27 )/*(/ 3 ) + 2As-^(i-^)/|(/3) + E*(/3, s) 


s^-y + 2\sy + b 

with & = u; + cVa(/9)- 

To invert the Laplace transform, we set 

1_5(2-27) 

(4.11) A(/3, s) = -^-- 

^ ^ ’ 5^7 2As'i -h b 


(4.12) 

(4.13) 


H(/3,s) = 
fi(^.s) = 


- 5 ( 1 - 7 ) 


-h 2Asi' -h b 
1 


2\sy + b 

Inverting the Laplace transform in fl4.1ip . we obtain 

s1-5(2-27) 


A{Ft) = C-^ 

1 


^-v 


§27 + 2As+ -h E 

c-^ 


( 

-gl-5(2-27) 

5-5(1-27)- 

1 


— r] 
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(4.14) = ^ [^7+2^-275-2 " E,,,+25-2,S-l{vt^))] , 

where r] and ^ are the roots of the quadratic equation 

+ 2XE + b = 0. 

The functions appearing in fl4.14p are the classical one parameter Mittag-Leffler 
functions (see for example [7] for a complete review). 

The other terms appearing in fl4.1Up can be handled in a similar way. In particular, 
we have that 

(4.15) SiP, t) = 
and 

(4.16) Sl(/3,f) = ^|B,.,(fC) - B,.,(r,(7]. 

Taking all these terms together and applying the convolution theorem, we have the 
Fourier transform of the solution of the Cauchy problem 

{E,,,^,S-2,5-im - i?,,,+2^-275-1 (h^"))] flW) 

+ [E-^ {E,^,+s-psm - E,,,+s-p5{vt^))] fm 
+ T'^~^E*{(3,t- r)[E^,j{^T'^) - .E.^,^(r7r^)]cir|. 

Finally by using the inverse Fourier transform we obtain the claimed result. □ 


Remark 4.2. Under the same conditions of Theorem f.l, if 6 = t), the space- 
fractional derivative becomes a Riesz derivative and the parameter b appearing in 
the solution assume the more simple form 

(4.17) h = u + E\(5\^. 

Moreover, for a; = 0 the solution obtained in this case reduces to the one given by 
Garg et al. in |5]. 


Remark 4.3. IFe observe that, for') = 1/2, 5 = 3/2, / 2 (x) = E{x) = uj = D, a = 2 
and fi{x) = d(x), the solution 04.71) coincides with 


(4.18) 


u{x, t) = 


“+00 


2cy/7ri J _ 
d , 


e 


—Alt; 


A 


A/n ( —VEw^ — I + 


dw 


A 


vh 


— X' 


d|x|<Clo} 


+ c [(5(x — cw) + 6{x + cw)] >dw, 


(see Theorem f.2 in [13] for the complete proof). This is the probability density of 
the telegraph process T{t), t > 0, with a Brownian time, that is 

W{t)=T{\B{t\)- 
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In the case 7 = 1, for all 6 G [0,1], f 2 ix) = F{x) = u = 0, a = 2 and fi{x) = d(x), 
the solution (I4TD coincides with the fundamental solution of the telegraph equation, 
that is the probability law of the telegraph process. Moreover, we observe that, in the 
case 7 = 1, (5 G [0,1], fi{x) = F{x) = oj = 0, a = 2 and f 2 {x) = S{x), the Fourier 
transform of the solution 04.71) becomes 

(4.19) u*{ft,t) = ^ - 

i-h 

which coincides with the characteristic function of the following distribution 

OO 

(4.20) P{T{t) G dx, IJ N{t) = 2k + l} = P{r{t) G dx, N{t) = 2k + 1}, 

k=0 

where N{t) is the homogeneous Poisson counting process, that governs the changes 
of direction of the telegraph process. 

Remark 4.4. The problem of the existence and uniqueness of the solution for the 
abstract fractional telegraph process has been considered by Herzallah and Baleanu 
m [8]. By simple adaptation of their arguments, the existence and uniqueness of the 
solution to the Cauchy problem in Theorem f.l can be proved. 
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